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The Walecka model contains essentially two parameters that are associated with the Lorentz 
scalar ( S ) and vector ( V ) interactions. These parameters are related to a two-body interaction 
consisting of S and V , imposing the condition that the two-body binding energy is fixed. We 
have obtained a set of different values for the nuclear matter binding energies ( Bn ) at equilibrium 
densities ( p 0 ). We investigated the existence of a linear correlation between Bn and p 0 , claimed 
to be universal for nonrelativistic systems and usually known as the Coester line, and found an 
approximate linear correlation only if V — S remains constant. It is shown that the relativistic 
content of the model, which is related to the strength of V — S , is responsible for the shift of the 
Coester line to the empirical region of nuclear matter saturation. 
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1. INTRODUCTION 

Quantum Chromodynamics (QCD) is the fundamental theory of the strong interaction and hence it 
should be able to explain possible modifications of hadron properties in the nuclear medium. However, 
typical nuclear phenomena at intermediate and low energies cannot be analitically derived from QCD, 
although one hopes that QCD will be solved numerically on the lattice in the near future. Meanwhile, we 
are left with the construction of phenomenological models in order to describe nuclear phenomena and 
bulk properties. On these grounds, we can conceive models in which hadrons are the degrees of freedom of 
some proposed Lagrangian density. From the Lagrangian, a two-nucleon interaction can be obtained and 
the parameters of this interaction are adjusted to reproduce the two-nucleon bound and scattering states 
observables. With this two-nucleon interaction, the many-body problem can be solved via a Brueckner- 
Hartree-Fock (BHF) calculation or a Relativistic-Brueckner-Hartree-Fock (RBHF) ASS , if one intends 
to incorporate relativistic effects. 

Another phenomenological way to treat the many-body problem has been proposed by Walecka and 
collaborators Hi Their approach is based on a simple renormalizable field-theoretical model which 
is often referred to as Quantum Hadrodynamics (QHD). In this model nucleons interact through the 
exchange of a and lo mesons, with the a providing medium range attraction and the u> the short-range 
repulsion. This model is usually solved in the mean field approximation (MFA), in which the meson 
fields are replaced by their expectation values. The Walecka model has achieved important goals in the 
description of hadronic matter such as, for example, the calculation of some bulk properties of nuclear 
matter as well as some properties of finite nuclei. An interesting result of the model was to show the 
relativistic mechanism for nuclear matter saturation: it occurs at a density ( p Q ) at which the scalar (S) and 
the vector (V) potentials largely cancel each other out. A curious aspect of this model is that the masses 
of the scalar and vector mesons ( m a and m u ), and the coupling constants ( g a and g w ) are eliminated 
in the equation of state for infinite nuclear matter in favour of C 2 = g 2 M 2 /m^ and C 2 = g^M 2 /in 2 , . 
These dimensionless constants, C 2 and C 2 , the only free parameters in this simple version, are fitted to 
reproduce the nuclear matter bulk properties. Besides the relativistic mean field calculation, the Walecka 
model has also been used in a more complete treatment, the relativistic Hartree-Fock approximation 
Sill A more sophisticated version of the Walecka model that includes non-linear self-interactions 
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of the scalar field has been proposed by Boguta and Bodmer @]. This non-linear version of the Walecka 
model has obtained not only more reasonable results for the incompressibility and the effective nucleon 
mass at nuclear matter saturation density but also has described various aspects of nuclear structure 
properties MM This model with a very limited number of parameters is also able to describe deformed 
nuclei m l| and for the first time the anomalous shifts in the isotopic chains of different nuclei has been 
explained jl3L 0- As a consequence of all this work, a new parametrization for this non-linear version 
of the Walecka model has been proposed and gives a very good description not only for the properties of 
stable nuclei but also for the nuclei far from the valley of beta stability E3- 

The RBHF and Walecka models treat the many-body problem in different perspectives. The RBHF 
approach starts from a two-nucleon interaction which, in principle, fits some of the two-nucleon ob¬ 
servables and proceeds iterating an effective two-nucleon t-matrix up to the point where convergence is 
obtained to the nuclear matter saturation. Therefore, depending on the two-nucleon interaction used, 
the nuclear matter binding energy and the saturation density change. On the other hand, the Walecka 
model and some of its variants try first of all to fit directly the observed nuclear matter binding energy 
at the experimental density. The tracing back of the two-nucleon interaction becomes difficult in this 
mean field scenario. To understand why, let’s point out that the two-nucleon interaction that emerges 
from the Walecka model in the static approximation is the sum of two Yukawa potentials, one attractive 
and the other repulsive. The range of these potentials should be controlled by the masses of the mesons, 
m a and m u . However, the mean field approach in a Hartree perspective shows a dependence only on the 
dimensionless constants C% and C as already pointed out, and we loose track of the value of m c , and 
m w . Let us recall that, in the mean field picture, this happens in a medium where each nucleon feels 
equally an averaged constant interaction coming from all the other nucleons. If one intends to have a 
crude idea about how a mean field many-body calculation such as the Walecka model may be connected 
to a two-body system in the vacuum, one has to model a two-body interaction in a relativistic perspec¬ 
tive. The simplest interaction should be the sum of two large potentials, a scalar one S and the time 
component of a vector one V , each of them with an absolute value of a few hundreds MeV. In infinite 
nuclear matter, in the framework of the Walecka model, these two components cancel each other out to a 
large extent in the calculation of the binding energy. The assumption is that the same mechanism applies 
in the two-nucleon system to calculate the two-nucleon binding energy. 

Calculations done for the nonrelativistic nuclear matter BHF model pj, performed with different two- 
nucleon interactions, all of them fitting the deuteron binding energy and scattering data, show a correla¬ 
tion between the binding energy and the corresponding saturation density. This correlation, known as the 
Coester line, shows that the results obtained using different nucleon-nucleon interactions predict nuclear 
matter saturation to be located along a line which does not include the empirical data [1611 . Recent work 
m based on a relativistic model of the nucleon-nucleon interaction addresses the discussion of whether 
this same correlation also takes place in a relativistic many-body calculation, exhibiting thus a kind of 
universality. We focus on the same question and, through a simple model, find that such a correlation 
is indeed present when the quantity V — S is approximately constant. This last assumption can be 
justified since results from different hadronic models suggest that, in order to reproduce the experimental 
spin-orbit splittings for finite nuclei, the V — S quantity should lie in a narrow band E3- Moreover we 
conclude that the relativity within the model, which is related to the strength of V — S , is responsible 
for the shift of the Coester line to the empirical region of nuclear matter saturation. We see this as an 
indication of a Dirac structure for the N-body potential with large scalar and vector parts. 

After having finished this work we became aware of the appearance of a very recent paper analyzing 
this same correlation with the Bonn potentials; these involve different types of meson exchanges, in a 
relativistic Brueckner approach 0. Representing the in-medium on-shell T matrix covariantly by five 
Lorentz invariant amplitudes, this last paper also arrives at the same conclusion as ours, namely that the 
new Coester lines are shifted towards the empirical region of saturation supporting the Dirac structure 
of the Nucleus-Nucleus potential. 
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Thus we can see that, despite dealing with a simple two-body model that has only central terms and 
models nuclear matter with the simple linear version of the Walecka model, we have the right physics, 
which is embeded in the relativistic form of the potential with large scalar and vector parts. 

2. THE TWO-BODY MODEL 

We start out with a two-body interaction consisting of a repulsive part ( V ) and an attractive one 
( S ). For simplicity, we will take the potentials as constant inside a determined radius R 0 and zero 
otherwise. The relativistic interaction is given by the operator U = Ur = S + 7 °V, where 7 0 is the 
time component of the Dirac matrix and S is implicitly multiplied by the identity matrix. We simulate 
our two-body system by a one-body Dirac equation in spinorial form, 

{Ej° — y-p — M — U) ip = 0. (1) 

The above equation can be solved for U = 7 °S ( the relativistic case where E is considered as the time 
component of a pure attractive vector interaction, hereafter labeled R1 ) or for U = Ur (the relativistic 
case with scalar and vector interactions, R2 ). We have solved Eq. m for each of these two potentials 
and for zero orbital angular momentum in coordinate space for different values of R 0 , varying from 1.35 
to 1.7 fm, obtaining different values for the two-body binding energy B -2 = E — M. Since we will later 
establish a connection with a N-body relativistic Walecka model, we have already imposed, for the case 
U = Ur , a constraint on S and V to satisfy the optical relationship required from the Hugenholtz-van 
Hove theorem for any relativistic nuclear matter model in the mean field approximation [20l . 

— B n = V + [kf + (M + S) 2 ] 1/2 - M, (2) 

where Bn is the nuclear matter binding energy per nucleon, kf the Fermi momentum at saturation 

density, and M the bare nucleon mass. We have identified in the relativistic model R2 the quantity 
5 + 17 with the depth of the square well vector potential E . This identification comes often when one 
interprets the Walecka model as qualitatively performing a large cancellation between S and V , resulting 
in a value close to —70 MeV for E. In Fig. (1) we display our results for the case B N = 16 MeV, 

kf = 1.3/m -1 , R 0 = 1.4 fm, and M= 939 MeV where, for the sake of comparison, a nonrelativistic 

Schroedinger equation calculation for the nonrelativistic interaction Unr = E is included. The curves 
show that the solution of the relativistic Dirac equation with the time component of a pure vector 
interaction changes the nonrelativistic Schroedinger result only slightly. The binding energy in the former 
has a value somewhat greater than in the latter. However, the Dirac equation solved with the relativistic 
interaction, containing scalar and vector parts, changes dramatically the two-body binding in the opposite 
direction, i.e. the system becomes less bound. 

The results can be understood as follows. The three distinct cases can be presented through formal 
dispersion relation equations in terms of the large component x of the wave function, 


p 2 x 

= 2 M (e - 

s)x, 

NR Model 

(3) 

p 2 x 

= 2 M ( e - 

S )( 1+e 2 

R1 Model 

(4) 

p 2 x 

= 2 M (e - 

S )< 1 +e 2 M >* 

R2 Model 

(5) 


where in the R2 model A = V — 5, E = 5 + V, and in Eqs. 0-0) e stands for the difference between 
relativistic energy and the nucleon rest mass. 

All these cases can be cast into an unified formula in terms of a nucleon effective mass M*, 



2 M 


X = (e - E)x, 


( 6 ) 
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where 


M* = 

M, 

NR Model 

(7) 

M* = 

M (1 + e ~ E ), 
v 2 M h 

R1 Model 

(8) 

M* = 

M < 1+C 2M >• 

R2 Model 

(9) 


Since £ and A are negative and positive quantities respectively, Eqs. ©-© suggest that the model 
R1 (where the effective mass is larger than M , leading to a smaller kinetic energy) should in fact lead 
to a more tightly bound system, with greater binding energy than model NR, whereas model R2 (where 
the effective mass is smaller than M , leading to a larger kinetic energy) is the case where the system 
is less bound. The same information could be inferred if one thinks in terms of the effective quantities 
£* = £(l-£/2M) and e* = e (1 — E/M) for Eq.©, and £* = £(l-A/2M)and e* = e(l-A/M) 
for Eq.©. In the cases we are discussing, the binding energy is small compared to the nucleon rest 
mass and we have neglected terms of order e 2 /2 M. Note that for the model R2, not only E but also 
A , appear in Eq. ©• The potential depth E and the radius R a control the two-body binding energy 
in Eqs. 0 and ©, whereas in Eq. © an additional value - A - is needed, working therefore as a 
new degree of freedom of the problem, since for the same R 0 , different possible values of E and A are 
possible, as shown in Fig. 2. The large difference between the results of the models R1 and R2 is not 
surprising, since A is expected to be about ten times larger than E. To stress this point, it is worth 
mentioning that the two-body bound state for the R1 model arises when — £f ? 2 > ir 2 ft 2 /4M while for 
the R2 model the condition reads — E A R 2 > ir 2 h 2 /8 , in units where c = 1 . 

3. THE N-BODY MODEL 

The N-body model we use for nuclear matter is the linear Walecka model that now we briefly present. 
The degrees of freedom are baryon fields (ip), scalar meson fields (a), and vector meson fields (qj). The 
Lagrangian density is given by 

C = ipi'y l _ l d^ip — ipMip + ^(d^ad^a — m 2 a 2 ) 

+ g s vipip - + im 2 ^^ - gy^^ip^, ( 10 ) 

where F fw = d^u , M is the bare nucleon mass, and m s and m v are the scalar and vector mesonic 
masses respectively. From the above Lagrangian we obtain, through the Euler-Lagrange formalism, the 
equations of motion for the nucleon and mesons fields. 

When the meson fields are replaced by the constant classical fields a 0 and u> 0 , we arrive at the mean-field 
approximation with the equations 


UJ 0 

<Jo 


9v 


9v 


;(lp + 1p) = p b , 

rriy rriy 

9 s / 7 / \ 9s 

— (iPiP) = —p s , 


mt 


m± 


( 11 ) 

( 12 ) 


where pb and p s are the vector and scalar densities respectively. 

Now we define the scalar ( S ) and the vector ( V) potentials. This can be done by looking at the Dirac 
equation for the models, and rewriting M* in the form M* = M — g s a = M + S. Still from the 
analysis of the Dirac equation, V can be defined as a quantity which shifts the energy, V = g v u 0 ■ It is 
convenient to introduce the following dimensionless quantities: m* = M* / M , C 2 = g 2 M 2 /m 2 , and 
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U 2 = glM 2 /m 2 v . Using Eqs. (EJ) and m we can rewrite the scalar and vector potentials in terms of 
Cl and Cl as 


S = 


cl 

M 2 Ps 


v _cl 

V M 2pb ' 


(13) 


The expressions for the energy density and pressure at zero temperature can be found as usual by the 
average of the energy-momentum tensor, 




2M 2 2 Cl 


27r 2 


k 2 dk E*(k), 


and 


a 


M 4 


ui 2 ly± /I *\2 | 1 7 f f I 2 

P= 1KP P + 32^l kdk wwr 


(14) 

(15) 


where p = ( 7 / 67 r 2 )/cj , 7 = 4 is the degeneracy factor for symmetric nuclear matter, and E*(k) = 
(k 2 + M * 2 ) 1 / 2 . 

The solution of the model is obtained explicitly through the minimization of £ as a function of m*. 
The equation thus obtained reads 


1 * 
1 — m 


7 r Xf x 2 dx 

2 tt 2 J 0 (x 2 +m* 2 )i 


(16) 


where we have introduced the dimensionless variable x = j-. This equation has to be solved self con¬ 
sistently and provides the basis for obtaining all kinds of thermodynamical quantities in the mean field 
approach we are using. 

The usual procedure to obtain the values of the coupling constants is the following: at the sat¬ 
uration density p 0 the pressure vanishes ( hydrostatic equilibrium ) and at this same point the 
Hugenholtz-van Hove theorem allows us to write down the relation, already built-in in Eq.(2), £/p a = 
„l 1/2 


V 


k 2 


(M + sy 


— M. These two constraints can then be used to find the constants C 2 and 


C 2 by imposing the experimental values for the nuclear matter binding energy £/p 0 — M = —Bn = -16 
MeV and for the density p 0 = 0.15 fm~ 3 . 

We could have used the non-linear version of the Walecka model discussed in the introduction 0 - 
This would have changed the expression for the energy density and the self consistent Eq. d; however, 
the correlation we are analyzing in this work would still survive, because of the presence of the large 
scalar and vector potentials. 


4. RESULTS AND DISCUSSIONS 

To look for the correlation between two-body and N-body bindings we proceed in the following way. 
From the two-body calculation for a fixed two-body binding energy B 2 at a given radius R 0 one obtains 
the scalar ( S ) and the vector ( V ) potentials. In the Walecka model we then use these two quantities, 
coming from the two-body calculation, in order to find the constants C 2 and C 2 . Only after we have 
these constants do we proceed in the calculation of the N-body system saturation. Therefore, we will 
have a set of saturated (Bn , p 0 ) pairs for each different ( S , V) set pairs, which keeps a fixed B 2 - We 
have chosen a fixed B 2 = 2 MeV, allowing variation of R 0 from 1.4 fm to 2.1 fm. The set of ( R a , S , V ) 
is given by Fig. 2 in a three-dimensional plot. 
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Within the above procedure we have obtained the set of (Bjy,p 0 ) values presented in Fig. 3. In 
principle, there is no clear correlation between Bn and p Q unless one constrains the interaction of the 
two-body system to have the same range R 0 but allowing different values of A, or to have the same 
value of A but having different R a values. In the first case, the correlation points show a positive slope 
for increasing saturation density while in the second one, the Coester line correlation is present. This 
correlation has been shown to exist in nonrelativistic calculations of the nuclear matter saturation point. 
Different phenomenological two-body potentials furnish different (Bn , p Q ) values, which lead to a rough 
linear increasing of Bn for increasing p a also shown in the Fig.3. To trace back the reason why this 
happens in the relativistic case, we have rewritten Eq. 13 as 


— B n 


J±_ 

2 M* 


(17) 


where M* = M (1 — ). Since Bn ■C A , M* —*■ M (1 — A/2 M) and when A becomes constant, 

Eq. 1171 1 acquires a nonrelativistic form with an effective nucleon mass. The relativistic information of 
the model is contained in the strenght of A . The relativistic effect is shown in Fig. 3, where the Coester 
line is shifted to the empirical region of saturation provided A is large. This suggests a Dirac structure 
for the potential with large scalar and vector parts. 

Regarding the role of A in the nuclear data fittings, it is interesting to point out that recent nuclear 
matter analysis with many different quantum-hadrodynamics models show a correlation between the 
finite nuclei spin-orbit energy splitting and the nuclear matter effective nucleon mass m* ll9l . The 
spin-orbit energy splitting increases as m* decreases. To accurately reproduce the empirical splittings, 
it is found that one should require m* to be between 0.58 and 0.64. From the Hugenholtz-van Hove 
theorem j2(|, we can relate m* and A through A = 2 M — Bn — M* — (k 2 + M* 2 ) 1 / 2 . Therefore, 
a correlation between the splittings and A also exists, namely the energy splitting increases with A. 
Moreover, the phenomenology of finite nuclei restricts the values of A to lie in a narrow band if one 
intends to reproduce good results for spin-orbit splitting. On the other hand, relativistic results, obtained 
via Dirac-Brueckner-Hartree-Fock calculations, usingthree different versions (A,B and C) of the Bonn 
potential, show a roughly linear Bn x p a correlation [l] . The corresponding nuclear matter A quantities 
lie on a narrow range. 

In Fig. 4 we present the incompressibility ( K) as a function of p 0 . Again, the possible correlations 
follow the same trend established in Fig. (3). Roughly, K increases linearly with p Q for fixed A values. 
When we consider bigger values for A , K increases, turning the equation of state stiffen In Fig. (5), the 
effective nucleon mass M* is displayed as a function of p 0 , also showing the same qualitative behavior 
for the possible correlations in terms of the two-body R a and A inputs. As expected by our previous 
discussion, when A increases M* decreases. In these two last figures we only want to exhibit the 
correlations between the incompressibility and the effective nucleon mass on one hand and the saturation 
density on the other. The specific values of the two first quantities are not important here and we would 
certainly get more reasonable results for them at the nuclear matter saturation density had we used the 
non-linear Walecka model H^j. 


5. CONCLUSIONS 

The plots of the binding energy versus the saturation density of nuclear matter that are obtained by 
nonrelativistic calculations with a large number of different two-nucleon interactions show a correlation 
that is known as the Coester line. In a nonrelativistic approach, the Coester line has to be conceived 
from the phenomenology contained in two-body systems with a low fixed binding energy interacting via a 
short range potential. In terms of our simple two-body model, the main ingredients are £ ( the depth of 
the potential) and the range R 0 . Both suffice to parametrize in a nonrelativistic approach the two-body 
problem itself, with a fixed two-body binding energy. 
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From a relativistic point of view, we have presented a discussion about this possible correlation. Our 
main conclusions can be summarized as follows: 

(1) In the relativistic case with scalar and vector interactions, a new degree of freedom A arises and 
we have seen that unless A is constant, the correlation contained in the nonrelativistic Coester line is 
destroyed. 

(2) Since the phenomenology of relativistic finite nuclei calculations restricts the values of A to lie 
in a narrow band, if one intends to reproduce good results for spin-orbit splitting, we conjecture that a 
Coester line would survive in such models. 

(3) For large values of A we show that the relativistic Coester line is shifted to the empirical region 
of saturation, which we see as an indication of a Dirac structure for the potential with large scalar and 
vector parts, and also as a signature of the relativistic content of nuclear matter. 

(4) For a fixed value of B 2 , the study of Bjv ,K and M* as a function of p 0 suggests that they are 
themselves correlated. 

We have modeled our system in a very simple way, where the two-body interaction has only the central 
term. We are aware that the phenomenological two-nucleon interactions in vacuum, as used by Reid, 
Paris, Hamada-Jonston, Bonn and others have not only simple central terms, but include spin-spin, 
spin-orbit and tensor force. Their magnitudes are different in strength. In particular, the tensor force is 
believed to be important for the nuclear matter saturation l2ltl . However, as already mentioned in our 
introduction, a recent paper has analyzed the Coester line for various Bonn potentials in a relativistic 
Brueckner approach [Ig]. Its conclusion is that the new Coester lines are shifted towards the empirical 
region of saturation, supporting a Dirac structure for the Nucleus-Nucleus potential. Then we can 
conclude that even when we are dealing with a simple two-body model that has only central terms, we 
have the dominant aspect of the correct physics. This physics is in fact embedded in the relativistic form 
of the potential, with scalar and vector parts, because, in the medium, many pieces of the interaction 
are averaged out, while the remaining part is largely dominated by the scalar and the vector central 
potentials psj . 

To finish, we remark that for the first time a connection between the Coster line has been proposed as 
a function of the A = V — S quantity, Fig. 3. We have exagerated the A interval in this figure to better 
exhibit how the relativistic content of the model would manifest itself in the Coester lines themselves. 
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FIG. 1: Two-body binding energy B 2 as a function of E = V + S for different approaches, Nonrelativistic (NR), 
Relativistic Dirac with pure vector interaction (R 1) and Relativistic Dirac with scalar and vector interaction (R'2). 



FIG. 2: Set of R = R 0 , S + V = E, and V — S = A values, with the constraint B 2 = 2 MeV , from the R2 model. 
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FIG. 3: N-body binding energy Bn as a function of the density po for different values of A. The curves correspond 
to A = 650 (A), 700 (B), 750 (C), 800 (D), and 850 (E) in MeV. The points on each line, from left to right, 
correspond to values of R 0 ranging from 1.4 to 2.1, in steps of 0.05 fm. All the points correspond to ( Ro, S , V ) 
values leading to B 2 = 2 MeV. The straight line represents the fitting of the nonrelativistic Coester line of ref. 1. 
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FIG. 4: Incompressibility K as function of the density p. A and R 0 are the same as in Fig. 3. 







